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Abstract. We establish a convenient necessary and sufficient con- 
dition for a multiaffine real polynomial to be stable, and use it to 
verify that the half-plane property holds for seven small matroids 
that resisted the efforts of Choe, Oxley Sokal, and Wagner [5]. 



In recent years, matroid theory has found connections with certain 
analytic properties of real multivariate polynomials. These properties 
are abstractions of physical characteristics of an electrical network. Not 
all matroids exhibit the same physically sensible behaviour that graphs 
do. It is an interesting (and often challenging) problem to determine 
whether a given matroid satisfies one or another of these physically- 
motivated conditions. In this paper we deduce a convenient necessary 
and sufficient criterion (Theorem 3(c)) for the "strong Rayleigh prop- 
erty" , and use it to verify this property for some small matroids, among 
them the Vamos cube V 8 . This supplements Branden's result [2] that 
the strong Rayleigh property is equivalent to the "half-plane property" , 
and resolves some questions left open by Choe, Oxley, Sokal, and Wag- 
ner [5]. 

Let Z(yi, ...,y m ) be a polynomial with real coefficients, and let E = 
{l,...,m}. The polynomial Z has the half-plane property (HPP) or is 
Hurwitz stable provided that whenever Re(y e ) > for all e G E then 
Z(yi, ...,y m ) 7^ 0. The polynomial Z is stable provided that whenever 
Im(y e ) > for all e G E then Z(yi, ...,y m ) ^ 0. Note that a homoge- 
neous polynomial Z is stable if and only if it is Hurwitz stable. If every 
variable y e for e G E occurs in Z to at most the first power, then the 
polynomial Z is multiaffine. For any index e G E let 

z. = ±z 

oy e 
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be the contraction of e in Z, and let 

Z e = Z\ n 

\y e =0 

be the deletion of e from Z. This notation is extended to multiple 
(distinct) sub- and super-scripts in the obvious way. For distinct indices 
e,f eE let 

AZ{e,f} = Z e Z f -Z ef Z. 

We refer to this as the Rayleigh difference of {e, /} in Z. If Z is 
multiaffine then Z = Z e + y e Z e for any e G E, and it follows that 

AZ{e,f} = Z{Z}-Z ef Z". 

The polynomial Z has the Rayleigh property provided that whenever 
y c > for all c E E, then for all {e, /} C E, AZ{e, /} > 0. The 
polynomial Z has the strong Rayleigh property provided that whenever 
y c eR for all ceE, then for all {e, /} C £, AZ{e, /} > 0. 

The motivation for these definitions is discussed in [5, 9]. The HPP 
and stability are studied in [1, 2, 3, 4, 5, 9], and the Rayleigh property is 
studied in [3, 6, 9, 10, 11]. Clearly the strong Rayleigh property implies 
the Rayleigh property. In [2], Branden proves that for multiaffine real 
polynomials the strong Rayleigh property is equivalent to stability. We 
give another criterion equivalent to these in Theorem 3(c). 

Now, let Z(yi, y m ) be a multiaffine real polynomial. To investigate 
the strong Rayleigh property, consider any pair {e, /} in E, and any 
third index g e E \ {e, /}. Since Z is multiaffine, AZ{e, /} is at most 
quadratic in y g . A short calculation shows that 

AZ{e, /} = Ay] + By g + C 

in which 

A = AZ g {e, /} = Zf g Z% - Z efg Zf 

and 

B = Zl°Z? fg + ZjZ{ g - Z?Zl f - Z efg Z^ 

and 

C = AZ 9 {e, /} = Z{ 9 Z) 9 - Z 9 ef Z ef9 . 

If Z has the strong Rayleigh property then for any real values of y c 
(ce£ \ {e, /, g}) the polynomial AZ{e, /} is nonnegative for all real 
values of y g . Since this quadratic polynomial in y g does not change 
sign, its discriminant is nonpositive: B 2 — 4AC < 0. In fact, this 
discriminant has a surprising feature that can be put to good use. 

Proposition 1. Let Z(yi, y m ) be a multiaffine polynomial, and let 
e, /, g be distinct indices. Then the discriminant of AZ{e, /} with re- 
spect to y g is symmetric under all permutations of the indices {e, /, g}. 
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Proof. With notation as in the previous paragraph, one calculates that 
the discriminant B 2 — AAC equals 

(zl 9 z%) 2 + (zfzQ 2 + (ZfZ° ef ) 2 + (Z efg Z^f 
-2Z{°Z) g ZjZ{ g - 2Zf°Z%ZfZ° f - 2Z e /Zf g Z?Z° ef 

-2(z£°z e fg + zfzl 9 + z? zl f )z efg z^ 

■\Zp'Z'/Z:/Z !f:l + AZ%Z[ g Z 9 ef Z^. 

This is clearly symmetric under all permutations of the indices {e, /, g}, 
as claimed. □ 

Lemma 2. Let Z(yi, ...,y m ) be a multiaffine polynomial with positive 
coefficients, and let e, f, g be distinct indices. Fix real values for all y c 
(ce E \ {e,f,g}). Assume that AZ f {e,g} > 0, that AZ f {e,g} > 0, 
and that AZ{e,f} > for all y g G R. Then AZ{e,g} > for all 

y f eR. 

Proof. We have AZ{e, g} = Ayj + Byj + C with real coefficients A > 
0, B, and C > 0. By Proposition 1, the discriminant of AZ{e,g} 
with respect to yj equals the discriminant of AZ{e, /} with respect to 
y g . By the hypothesis of the Lemma, this discriminant is nonpositive. 
Consequently, AZ{e,g} does not change sign for yf e R. It follows 
that if either A = or C = then B = 0. Therefore, AZ{e, g} > for 
all yj £ R, as was to be shown. □ 

Theorem 3. Let Z(yi, ...,y m ) be a multiaffine polynomial with positive 
coefficients. The following conditions are equivalent: 

(a) Z is stable; 

(b) Z has the strong Rayleigh property; 

(c) for every index e G E, both Z e and Z e have the strong Rayleigh 
property, and for some pair of indices {e, /} C E, AZ{e, /} > for 
ally c eR (ceE\ {e,/}). 

Proof. The equivalence of conditions (a) and (b) is Theorem 5.6 of 
Branden [2]. That (b) implies (c) is immediate, since the operations 
Z i— > Z e and Z i— > Z e both preserve the strong Rayleigh property (see 
Proposition 5.4 of [6]). It remains only to show that (c) implies (b). 
Assume the hypothesis of (c). Then AZ{e,f} > for all y c G R 
(c G E \ {e, /}), by hypothesis. By one application of Lemma 2, we 
deduce that for any g G E \ {e,f}, AZ{e,g} > for all y c £ R 
(c G E \ {e, #}). The cases for AZ{f,g} follow by symmetry. For 
a pair {g, h} disjoint from {e, /} we apply Lemma 2 again, beginning 
with AZ{e,g} > for all y c G R (c G £ \ {e,#}), to deduce that 
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AZ{g,h} > for all y c G R (c G E \ Thus, Z is strongly 

Rayleigh. □ 

The advantage of the criterion (c) is that only one Rayleigh difference 
needs to be checked, instead of the (™) that are required a priori. The 
strong Rayleigh condition on the minors Z e and Z e can be assumed if 
one is proceeding by induction on \E\. We give a few applications of 
this method below. 

The above concepts apply to a matroid M with ground-set E through 
its basis- generating polynomial 

M(y) = $> B , 

B 

in which the sum is over the set of all bases B of M, and y B = Yl e &B He- 
These polynomials are clearly homogeneous and multiaffine. Thus, 
M(y) has the HPP if and only if it is strongly Rayleigh (Corollary 
5.14 of [2]). Notice that M e is the basis-generating polynomial of the 
deletion M \ e, and that M e is the basis-generating polynomial of the 
contraction M/e. The Potts model partition function of M is a much 
more informative polynomial to which some of the above conditions 
can also be applied - see [8, 11]. 

The study of matroids with the HPP was begun in [5], from which 
we require the following results. 

• The class of HPP matroids is closed by taking duals and minors. 

• Every matroid with at most six elements is HPP. 

• Every matroid with rank (or corank) at most two is HPP. 

• Up to duality, there are only four matroids with seven elements for 
which the HPP is undetermined: these are ^f 4 , W 3+ , W 3 + e, and J" 7 
depicted in Figure 1. They all have rank three. 

• The rank three seven-element matroids 5Ff 5 and 7" are HPP. (The 
precise structure of these matroids is not important here - see [5].) 

The HPP remained unsettled in [5] for three more small matroids. 
These are the Vamos cube V 8 (of rank 4) and the two one-point dele- 
tions nty \ 1 and nD 3 \ 9 (of rank 3) of the non-Pappus matroid n7 
depicted in Figure 2. (The non-Pappus matroid itself is not HPP.) 
Using the criterion of Theorem 3(c) we now show that all seven of 
these matroids do have the HPP/strong Rayleigh property. To show 
that AM{e, /} > for all y c G R (c G E \ {e, /}), we write AM{e, /} 
as a sum of squares of polynomials. (Though not needed for the present 
purposes, when doing this one should be aware of the theory developed 
around Hilbert's 17th problem; see [7], for example.) 
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3 4 5 3 4 5 3 4 5 3 4 5 



'F- 4 W 3+ W 3 +e <TL 

/ / 

Figure 1. The four undetermined cases with \E\ < 7. 




Vo nT 



Figure 2. The Vamos cube and the non-Pappus matroid. 

Every one-point deletion or contraction has six elements, so is strongly 
Rayleigh by the above remarks. A short calculation verifies that 

AM{1,2} 

= ^(2/32/7 + 2/52/7 + 2/42/5 + 2/32/4 + 2/32/5 + 2/32/6 + 2/62/7 + 2/42/e) 2 
+ ^(2/32/7 + 2/32/4 + 2/32/6 + 2/32/5) 2 
+ ^(2/42/6 - 2/S2/7) 2 + \{y&Vi ~ 2/42/s) 2 
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By Theorem 3, it follows that 3^ is strongly Rayleigh. 

• M = W 3 +: 

Every one-point deletion or contraction has six elements, so is strongly 
Rayleigh by the above remarks. A short calculation verifies that 

AM{1,2} 

= ^ (2/62/3 + 2/42/7 + 2/32/5 + 2/52/7 + 2/42/6 + 2/62/7 + 2/42/3 + 2/42/5) 2 

+ ^ (2/62/3 + 2/32/5 + 2/42/3 + 2/42/5 ) 2 + ^ (2/42/7 + 2/52/7 + 2/62/7) 2 + ^2/42/6 
By Theorem 3, it follows that W 3+ is strongly Rayleigh. 

• M = W 3 + e: 

Every one-point deletion or contraction has six elements, so is strongly 
Rayleigh by the above remarks. A short calculation verifies that 

AM{1,2} 

= g + y3y6 + y4V7 + ^ 6 + y3V4 + yey7 + yhyi + y3V5 + ^2/s) 2 

+ ^(2/32/7 + 2/32/6 + 2/32/4 + 2/42/5 + 2/32/5) 2 + ^ (2/42/6 - 2/52/7) 2 

1 1 

1 2 2 1 - 1 - 2 2 

+ 22/42/7 + 2^6l/7 

By Theorem 3, it follows that W 3 + e is strongly Rayleigh. 

• M = r 7 : 

Every one-point deletion or contraction has six elements, so is strongly 
Rayleigh by the above remarks. A short calculation verifies that 

AM{1,2} 

= 2 (y*vt + y4ye + 2/42/7 + V3Va + y5y7 + y5m + y6m + yiy ^ 2 

+ ^(2/32/4 + 2/32/7 + 2/52/3 + 2/62/3 + 2/42/5) 2 

+ ^(2/42/6 -2/S2/7) 2 + 2-2/42/7- 
By Theorem 3, it follows that 7' 7 is strongly Rayleigh. 

• M = n7 \ 1: 

Every one-point contraction has rank 2, so is strongly Rayleigh. Every 
one-point deletion is isomorphic to one of 4 , W 3 + e, ^Ff 5 , 3" 7 , or 
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7". By the results above, these are all strongly Rayleigh. A short 
calculation verifies that 



AM{2,4} 

= ^ (2/32/7 + 2/32/9 + 2/52/9 + 2/62/7 + 2/72/8 + 2/72/9 

+2/52/7 + 2/52/8 + 2/62/8 + 2/82/9 + 2/32/6 + 2/32/5) 2 

+ 2 (2/62/7 + 2/72/8 + 2/72/9 + 2/32/7 + 2/52/7) 2 

1 11 

+ 2 (2/32/9 + 2/52/9 + 2/32/5 - 2/62/8 ) 2 + 2(2/52/8 - 2/32/6 ) 2 + 2^2/9 



By Theorem 3, it follows that n7 \ 1 is strongly Rayleigh. 
• M = n?\9: 

Every one-point contraction has rank 2, so is strongly Rayleigh. Every 
one-point deletion is isomorphic to one of S^ 4 or 7' 7 . By results above, 
these are both strongly Rayleigh. A short calculation verifies that 



AM{1,2} 

= ^ (^ 6 + + y z yi + yeV7 + ^ 8 + 2/52/s 

+2/32/6 + 2/32/5 + 2/42/5 + 2/52/6 + 2/42/8 + 2/32/4) 2 
+ 2 (2/32/8 + 2/32/6 + 2/32/7 + 2/32/5 + 2/32/4 + 2/42/8 ) 2 

+ 7^(2/42/5 + 2/42/6 + 2/52/6 - 2/72/s) 2 + 7^(2/52/8 - 2/62/7) 2 



By Theorem 3, it follows that nT \ 9 is strongly Rayleigh. 
• M = V 8 : 

The Vamos cube is self-dual. Every one-point contraction is isomorphic 
to ^Ff 4 or jFf 5 . Every one-point deletion is isomorphic to (3 r 7 ~ A )* or 
(jFf 5 )*. By the results above, these are all strongly Rayleigh. A short 
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calculation verifies that 
AM{1,2} 

= \ (2/32/42/5 + 2/32/42/6 + 2/32/42/7 + 2/32/42/8 + 2/32/52/7 + 2/32/52/8 
+2/32/62/7 + 2/32/62/8 + 2/42/52/7 + 2/42/52/8 + 2/42/62/7 

+2/42/62/8 + 2/52/62/7 + 2/52/62/8 + 2/52/72/8 + 2/62/72/s) 2 

+ 1 (2/32/42/5 + 2/32/42/6 + 2/32/42/7 + 2/32/42/8 

+2/32/52/7 + 2/32/52/8 + 2/42/62/7 + 2/42/62/s) 2 

+ 1 (2/32/42/5 + 2/32/42/6 + 2/32/42/7 + 2/32/42/8 

+2/32/62/7 + 2/32/62/8 + 2/42/52/8 + 2/42/52/7) 2 

+ 1 (2/32/42/5 + 2/32/42/6 + 2/32/42/7 + 2/32/42/8 

-2/62/72/8 - 2/52/72/8 - 2/52/62/7 - 2/52/62/s) 2 



+ ^(2/32/62/7 


- 2/32/52/8 


+ 


2/42/62/7 




2/42/52/8 


+ 


2/62/72/8 - 


2/52/62/s) 2 


+ ^(2/32/52/7 


- 2/32/62/8 




2/42/62/8 


+ 


2/42/52/7 


+ 


2/52/72/8 - 


2/52/62/s) 2 


+ ^(2/32/52/7 


+ 2/32/62/7 


+ 


2/42/62/8 


+ 


2/42/52/8 


+ 


2/52/62/7 + 


2/52/72/s) 2 


+ ^(2/32/52/8 


+ 2/32/62/8 


+ 


2/42/62/7 


+ 


2/42/52/7 


+ 


2/52/62/7 + 


2/62/72/s) 2 


+ ^(2/32/62/7 


- 2/32/52/8 


+ 


2/42/62/7 




2/42/52/8 


+ 


2/52/62/7 - 


2/52/72/s) 2 


+ ^(2/32/52/7 


- 2/32/62/8 




2/42/62/8 


+ 


2/42/52/7 


+ 


2/52/62/7 - 


2/62/72/s) 2 


+ ^(2/32/52/7 


+ 2/32/62/7 


+ 


2/42/62/8 


+ 


2/42/52/8 


+ 


2/52/62/8 + 


2/62/72/s) 2 


+ ^(2/32/52/8 


+ 2/32/62/8 


+ 


2/42/62/7 


+ 


2/42/52/7 


+ 


2/52/62/8 + 


2/52/72/s) 2 



By Theorem 3, it follows that V 8 is strongly Rayleigh. 
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